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ABSTRACT 


The  solution  that  is  developed  herein  can  be  used  to  compute  thermal 
stresses  and  displacements  in  the  moderately  thick  walls  of  appreciably 
curved  ogival  radome  shells  whose  material  properties  change  with  heat 
intensity  along  the  span  and  through  the  wall.  No  restriction  is  placed 
on  temperature  dependent  variations  of  the  thermal  strain  expression. 

A  restriction  imposed  on  the  theory  is  that  the  wall  material  must 
not  be  overstrained.  Overstrain  is  any  amoun  beyond  the  elastic  limit 
in  excess  of  allowable  strain  that  accompanies  a  maximum  stress  equal  to 
the  usable  strength  of  the  material  as  determined  for  the  worst  heating 
condition  in  a  given  application. 

Three  examples  of  the  computerized  results  are  presented.  The  first 
example  deals  with  a  28.3"  x  6.75"  bicentrlc-oglve  radome;  the  second, 
with  a  thick-walled  hemispherical  radome;  and  the  third,  with  a  thin-walled 
hemisphere.  In  these  calculations,  linear  spanwlse  temperature  distributions 
are  v.sed  that  were  obtained  from  experimental  and  heat-transfer  analyses 
of  a  radome  tested  in  a  wind  tunnel.  Through  the  wall,  over  specified 
layers,  the  thermal-strain  function  is  approximated  as  being  linear  with 
temperature.  However,  the  temperature  varies  as  a  third-order  nonlinear 
function  through  the  thickness  of  the  wall  as  found  for  the  test  radome. 

It  appears  that  errors  due  to  the  approximation  amount  to  less  than  four 
percent  in  general. 

It  was  found  that  stresses  in  the  bicentric-ogive  radome  were  within 
safe  limits  in  regard  to  strengths  of  the  wall  material.  Good  agreement 
was  obtained  between  computed  and  experimentally  determined  stresses  at 
the  one  point,  1.75  inches  behind  the  tip,  where  strains  were  measured. 
Analyzed  for  comparison,  the  thick-walled  hemispherical  radome  was  over¬ 
stressed,  but  the  thin-walled  hemisphere  was  sufficiently  strong. 

It  is  also  concluded  that  additional  information  is  needed  on  tem¬ 
perature  distributions,  nosecap  behavior,  thermal  stress  maxima,  and  the 
effects  of  radome -missile  joints. 
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I .  INTRODUCTION 


In  his  fam  )us  treatise  on  elasticity,  A.  E.  H.  Love  [l]  discussed  the 
analysis  of  appreciably  curved  and  moderately  thick-walled  shells.  He  then 
proceeded  with  examples  In  which  he  made  certain  approximations  that  let  him 
neglect  terms  found  to  be  unimportant  within  the  range  of  consideration.  Many 
authors  thereafter  adopted  Love's  approximations  and  developed  the  modern 
thin-shell  theory,  [2],  [3],  [A]  ,  and  [s],  for  example. 

Today's  requirements  for  heat-resistant  structures  of  comparatively  thick- 
walled  shell  construction  have  led  to  further  investigation  of  the  problem. 
McDowell  and  Sternberg  [6]  examined  spherical  shells,  and  several  wr  t.’'s 
studii  d  cylindrical  shi  11s:  Tsao  [?],  Lee  [lO],  Hoff  and  Madsen  [li_, 

Rivello  [12]  and  [I6],  and  Bljlaard,  Dohrmann  and  Duke  [I8].  Conical  shells 
were  analyzed  by  Lu  and  Chang  [l4]  and  Weiss  [l7j.  Buckling  effects,  more¬ 
over,  were  included  in  [ll]  and  [1A],  and  sandwich  constructions  were 
discussed  in  [l2]  and  [I6]. 

Thermal  stresses  in  beams  were  computed  by  Barrekette  [S],  and  radome 
materials  for  high  temperature  usage  were  compared  by  Weckesser,  Hallendorff 
and  Suess  [is].  The  stiffness-matrix  method  was  employed  by  Dailey  [l3]  to 
obtain  numerical  solutions. 

In  [9],  F.  Lane  derived  solutions  for  shells  of  revolution,  and  D,  E. 
Magnus  and  D.  Eisen  presented  digital-computer  programs  for  numerical 
evaluations.  Lane  used  variational  calculus  to  minimize  the  strain  energy 
of  a  shell  to  obtain  his  results. 

Numerical  and  theoretical  analyses  developed  in  the  ensuing  text 
require  that  the  radome  wall  material  must  not  be  overstrained.  Overstrain 
refers  to  a  strain  that  exceeds  the  elastic  limit  at  a  concurrent  temperature. 
Reports  [l9]  to  [27],  inclusive,  were  prepared  as  part  of  the  overall 
thermal-stress  investigation. 


Numerals  in  brackets  denote  references  listed  at  the  end  of  this  report. 
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III.  NOMENCUTURE 


A,B,C 

D 

E 

F 

H 

I 

J,K 

L 

M 

N 

P 

Q 

R 

S 

T 

V 
W 

Y 
Z 


Constants  of  integration 

Wall  stiffness  factor,  D  =  Eh® /12(l-v^) 

Young's  modulus  of  elasticity,  psi 

Force;  Function;  Fahrenheit 

Intercept  height,  inches 

Wall  moment  of  inertia,  I  =  h®  /12 

Constants 

An  operator 

Moment  resultant,  ippi 

Normal-stress  resultant,  ppi 

Function 

Shear-stress  resultant,  ppi 
Spatial  radius,  inches 
A  series  or  its  sum 
Temperature,  °F 
Wall-slope  change,  radians 
Function 
Variable 

Geometric  axis  of  radome 


a  ,b 
c 
e 
f 
h 
i 
k 
I 

m,n 

P 

q 

r 

s 

u,w 

X 

y 

z 


Series  coefficients 

Half  thickness  of  wall  (c=h/2).  Inches 

Normal  strain,  ipi 

Function 

Wall  thickness,  inches 
Positive  integer;  1,2,3,... 

Constant 

Length,  inches 

Positive  integers;  1,2,3,... 

Function;  Pressure,  psi 
Auxiliary  shear  function,  ppi 
Planar  radius,  inches 
A  series  or  its  sum 
Displacement  functions,  inches 
Thickness  variable  ratio,  x  =  y/c 
Thickness  coordinate,  inches 
Coordinate  along  geometric  axis,  inches 


e 

e 

e 

V 

a 

T 

X 


Wall  bending  parameter,  radians 
Normal  strain,  ipi 

Coordinate  angle  of  rotation,  degrees  or  radians 
Poisson's  ratio 
Normal  stress,  psi 
Shear  stress,  psi 

Curvature  change,  radians  per  inch 
Coordinate  angle  of  azimuth,  degrees  or  radians 


The  following  symbols  are  used  as  subscripts; 


a:  Anterior  or  outer  surface  of  wall 

b:  Base  of  radome  at  the  body  of  a  missile 

c:  Central  surface  of  wall 

g:  General 
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1  ,m,n : 
o : 
p: 

q: 

B : 
t : 

y: 

r,0,z: 

R,e,\|i: 


Indices 

Origin  or  Initial  value  (zero) 
Particular 

Referred  to  function  q 
Secondary  or  inner  surface  of  wall 
Thermal  or  temperature 
In  the  y  direction 
Cylindrical-coordinate  directions 
Spherical-coordinate  directions 


IV.  THERMAL  STRESS  EXPRESSIONS 


1.  RADOME  WALL  STRESSES 


When  the  radome  mainbody  is  an  ogival  segment,  such  as  the  one  Illustrated 
in  Figure  la,  meridian  (a^)  and  hoop  (a@)  stresses  can  be  determined  with  a 
computer  program  based  on  the  following  equations. 


a 


a 


t’ 


(1) 


In  the  foregoing  expressions,  represents  a  function  having  the  dimenf  ons 
of  stress  (psl)  given  by 


^t  “  E€^/(1-v) 


(2) 


and  cjg  and  were  calculated  from  the  equations  written  below. 


R^a^  vyHCag-VCTj) 

~R  ■  (l-v2)rR 


a  = 


yH(ag-va^) 
R  "  (l-v2)rR 


(3) 


M,y 


=  TT  +  ~T~  > 
^  h  I 


M^y 


(4) 


Dimensions  h,  H,  R^.  and  variables  y,  r,  R  are  defined  on  Figure  la,  which  shows 
the  wall  profile  of  the  radome.  Symbol  "I"  denotes  moment  of  inertia  per  unit 
length  of  a  section  of  the  wall  (I  =  hP/12),  And  v  is  Poisson's  ratio  while 
normal  stress  (Nj^)  and  moment  (M^)  resultants  are  described  in  the  next  article. 
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2.  STRESS  AND  MOMENT  RESULTANTS 

Two  normal  stress  resultants  (I^  and  )  appear  in  equations  (4).  By 
definition,  they  are 

+  <5) 

where  all  four  of  the  right-hand  terms  are  obtained  from  the  following  integrals. 


Furthermore,  the  shearing  stress  resultant  (Figures  lb  and  2) 

+c 


-c 


(8) 


is  closely  related  to  the  first  two,  equations  (6),  as  discussed  in  Appendices 
A,  C,  D,  G,  and  References  [24],  [26],  and  [27]. 

Moments  that  bend  the  wall  of  the  radome  occur  in  equations  (4)  as  Mj^  and 
.  They  are  defined  to  be 

where  the  four  right-hand  terms  are  calculated  with  integrals  that  are  written 
below.  (Also,  see  Figure  lb). 
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-c  -c 


The  preceding  functions  are  evaluated  from  relationships  that  are  presented  in 
the  appendices  and  in  References  [l9]  to  [27],  inclusive.  Outer  8urfac_  tem¬ 
peratures  are  shown  in  Figure  3. 

3 .  SOLUTIONS 

(a)  General  Solutions  -  For  numerical  work,  it  is  frequently  expedient  to 
restrict  the  expansion  interval  of  infinite  series  that  are  employed  in  general 
solutions.  This  was  done  in  the  example  of  Part  VI-1  wherein  three  intervals 
span  the  length  of  the  radome.  The  nosecap  lies  within  0  ^  ^  »  64“  19  *23": 

(Figure  2),  and  the  mainbody  (Figure  1)  is  composed  of  two  segments  (MBl  and 
MB2)  whose  defining  intervals  are 


64'’19'23"  =  ijtjS  ^  ,|^=  73'’44'23"  (for  MBl) 

73°44'23"  =  \1(3  ^  ^  =  90°  (for  MB2) 

where  90°  coincides  with  the  radome  base  at  which  point  it  is  attached  to  the 
cylindrical  body  of  a  missile. 

For  general  solutions,  infinite  series  are  of  the  forms 


s 

m 


S 

ra 


00 


na  Y 
mn 


n-1 


n=o 


(12) 


with  numerical  results  discussed  in  the  appendices,  and  the  variable  Y  is  shown 
below  for  MBl  and  MB2. 
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Y  = 


sin  i|i  -  sin 
sin  ijig-  sin  * 

sin  iji  -  sin  \|(^ 
sin  i|(^  -  sin  * 


^  ^  +3 


^  i  ^  S  90‘ 


(13) 


(14) 


From  both  expressions  (13)  and  (14),  values  are  within  -1  ^  Y  ^  +1  from  left 
to  right  ends  of  each  interval.  Subsequently  (Appendix  D) ,  it  is  demonstrated 
that  the  general  solution  is 


N,  =  cos  it 

♦g 


m  m  ^ 


r 

Nfl  =  N,  +  ^ 

eg  ijtg 


m=l 


4 

z 

m*l 


mm 


(15) 


where  the  8^^^  (m  =  1,2, 3, 4)  are  constants  of  integration  determined  by  boundary 
conditions  as  explained  in  Appendix  E. 


The  functions  f^,  Fjj,  depend  on  Sj„,  being  computed  with 


f  (it)  1  „  *  s  ,  f  (*)  »  ,  ■  8  cos  t|t 

m  m=l,2  m  m  ’  m“3,4  m 


(16) 


F  (tjt)  10  =  k.S  cos  ijt,  F  (if)  .  -  k  S  cos^l  -  S  sin  t|t  (17) 
m  ^  m»l,2  i  m  m'^  m-3,4  i  m  ^  m 

wherein  represents  k^  in  MBl  and  kg  in  MB2  as  shown  below  in  equations  (18) 


k.  *  - - -  k,=  - 

^  sin  tjfg-  sin  ^  sin  \|t^  -  sin 


(18) 


The  bending  moments  are  calculated  from  changes  in  curvature  of  the  wall 
of  the  radome  with  orthogonal  components  expressed  by 


_ ^  ^  Y  -  ^  *^05  it 

^it  "  R  dit  ’  "  r 

^  c  c 


(19) 


where  V  denotes  change  of  slope  of  the  radome  wall  induced  by  thermoelastic 
effects.  The  moments  are 
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Mi=  -D(x^  +  vXg)v  Ms=  +  Xq) 


(20) 


D  =  EI/(l-v^)  =  Eh  / 12(1- 


(21) 


and  the  curvature  changes  are  calculated  from  the  next  two  equations. 


^0g 


_  cos  \lf 


y 


m=l 


mm 


’‘♦g  '  ’‘eg  ^  F 


m  m  ^ 


(22) 


m=l 


The  preceding  functions,  constants 
coefficients  a^^j^  ir.  equations  (12) 
described  in  great<.'r  detail  in  the 


and  and  computation  of  series 
by  means  of  recurrence  formulas  are 
appendices  and  in  Reference  [26], 


(b)  Particular  Solutions  -  In  general,  particular  solutions  cannot 
be  obtained  in  clcised  form  to  compute  stress  and  moment  resultants  in  the 
mainbody  (Figure  1).  It  is  demonstrated  in  Appendix  D,  however,  that 
approximate  formulas  can  be  developed  that  describe  the  desired  functions 
within  about  one- half  percent.  The  approximation  functions  are 


rBs  RjB^+BgCost) 


r  tan  A 
c 


Q  = 

.  r 


qR^sin  f 


,  =  Q^cot  ,!,  (23) 


N 


©p 


f  R  Bb 

-R.- 


cos  vli  -  — 
^  r 


K  cos  iji  ] 

csc^i  H - - -  (S,  +  B,  cos  \li) 

^  r  tan  ^  ^  ^ 

c  '  I 


(24) 


^  =  n- -R 


r  R  (B;+  B^cos  ^f) 


V  cos  \!i 


/ 

!  r- 


r  tan  \lr 
c 


^6  "  r 


(25) 


r  '  R 
c  c 


L 


R  cos  ill 

-  2  C  ^ 

cos  ill  +  —  CSC  'll  H - : - r 

’  r  tan  it 
c  ^ 


(B.J,'+  Bg'  cos  'll)  (26) 


in  which  q  /is  an  auxiliary  shear  function  used  to  find  the  total  shear  (Q.) 
by  means  of  equation  (23).  Integration  constants  B^jj  and  B^^  are  calculated 
in  Appendix/.  D  from  temperature  relationships,  and  they  are  fully  developed 
in  Referen«:e  [24]. 
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V.  THERMOELASTIC  DISPLACEMENTS 


1.  WALL  EXTENSION 


In  Appendix  F,  the  wall-extension  component  (w)  of  displacement  is 
developed  as 


w  =  w.sin  i  +  w  -  w^ 
i  q  t 


(27) 


where  w^^  is  a  constant  of  integration  (equalling  in  MBl  and  Wg  in  MB2), 

Wg  is  a  function  that  depends  on  the  auxiliary  shear  function  (q),  equations 
(ibO  and  (101),  and  Wj.  depends  on  the  temperature  distribution  as  given  by 
equations  (29),  (30),  and  (31). 


2. _ RADIAL  DISPLACEMENT 


After  finding  w  from  equation  (27),  the  radial  component  of  displacement 
(u)  is  calculated  with 


u  =  r^csc  i|i(l^  -  )/Eh  -  w  cot  i(t 


(28) 


where  the  coordinates  are  shown  in  Figures  la  and  lb,  and  the  other  terms  are 
as  defined  previously.  This  relation  for  u  is  obtained  in  Appendix  F. 


VI.  EXPERIMENTAL  DATA  AND  HEAT  TRANSFER  ANALYSES 
1 .  TEMPERATURE  DATA 

(a)  Wall-Thickness  Distributions  -  Available  information  on  radome  tem¬ 
perature  distributions  were  analyzed  and  reported  in  Reference  [22].  Tem¬ 
perature  distributions  that  are  shown  in  Figure  4  were  developed  in  the 
referenced  investigation  together  with  References  [20]  and  [2l]. 

(b)  Spanwise  Distributions  -  The  spanwise  temperature  distributions  found 
in  Reference  [22]  along  the  outer,  central,  and  inner  surfaces  of  the  test 
radome  are  plotted  in  Figure  5.  They  are  based  on  the  general  relationships 
developed  in  Reference  [20],  and  the  coefficients  are  empirical. 

(c)  Temperature-Distribution  Function  -  With  but  slight  error,  it  was 
found  in  Reference  L22J  that  the  approximate  temperature  distribution  through¬ 
out  the  wall  of  the  radome  was 


T  =  T  f(x) 
a 


(29) 
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DISTANCE  FROM  TIP  (lnehe») 

Fig.  5  EXTERIOR.  CENTRAL,  AND  INTERIOR  SURFACE  SPANWISE  TEMPERATURE 
DISTRIBUTIONS  FOR  28.3-INCH  x  6.7WNCH  VON  KARMAN  TEST  RADOME 
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where  T  “  T(\lf,y)  is  a  separable  function  of  and  x  =  y/c,  with  y  being  the 
wall  thickness  variable  and  \(i  the  coordinate  angle  as  in  Figure  la.  For 
thickness  variation,  we  have 

3 

f(x)  =  f^x"  =  0.257  +  0.358  x  +  0.305  x^  +  0.084  x^  (30) 

n=o 

in  which  numerical  values  of  coefficients  (f^^)  are  stated  in  the  expanded 
form  on  the  right. 

The  exterior-surface  temperature  distribution  is  given  by 


=  T^(t)  (31) 

which  is  shown  in  Figure  3  for  the  bicentric-ogive  radome  of  the  first  example 
in  Part  VIII. 

2 .  STRENGTH  DATA 

(a)  Strain  Measurements  -  Strains  were  measured  along  the  inner  surface 
of  the  wall  at  1.75  inches  from  the  tip  of  the  radome  emplc;ed  in  the  wind 
tunnel  tests  reported  in  References  [21]  and  [22].  No  breakage  was  reported 
for  test  results  used  in  analysis. 

(b)  Stress  Results  -  Stresses  were  calculated  from  the  strain  measurements. 
As  a  point  of  comparison,  an  experimentally  determined  meridian  stress  of 
18,000  psi  was  found  at  the  inner  wall  surface  at  1.75  inches  from  the  tip 

of  the  test  radome. 


VII.  NOSECAP  AS  A  SPHERICAL  SEGMENT 


1.  WALL  STRESSES 

Frequently,  a  radome  nosecap  is  the  spherical  segment  of  a  half-caliber 
ogive  pictured  by  Figure  2.  In  such  cases,  formulas  (3)  reduce  to 


a  =  R  a  /R, 
a  o  1  ’ 


(32) 


and  the  other  relations  are  unchanged.  Radius  Rq  is  the  principal  radius  of 
curvature  of  the  central  surface,  which  is  shown  in  Figure  2  also. 
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2 .  STRESS  AND  MOMENT  RESULTANTS 


The  normal  stress  resultants  (I^  and  Ng)  that  appear  in  equations  (4)  are 
given  by 


Ni=N^+N^.  N,=  N0+N^ 

where  the  right-hand  terms  are  obtained  from  the  following  integrals. 


(33) 


N 


N 


t 


(34) 


(35) 


Also,  the  wall  bending  moments  (1^  and  )  that  occur  in  equations  (4) 
are  defined  to  be 


(36) 


where  the  right-hand  terms  are  computed  with  the  integrals  shown  below. 


M 


+c 


R 

o 


ydy, 


-c 


(37) 


M 


t 


(38) 
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The  preceding  functions  are  evaluated  from  relationships  that  are  presented 
in  the  appendices,  and  outer  surface  temperatures  are  shown  in  Figure  3. 

3.  SOLUTIONS 

(a)  General  Solutions  -  For  the  nosecap,  whose  interval  is  shown  in 
Figure  2  as 


0  ^  ijt  ^  il(j=  64“19'23" 


the  general  solutions  corresponding  to  series  (12)  are  obtained  with  variable 
Y  replaced  by  1-cos  \|t ;  e.g., 


s 

m 


09 

Zb  (1-cos  i)”, 
mn 

n=o 


S  =  \  nb  (1-cos  ili)^ 

m  /  mn  ^ 

n=o 


(39) 


and  the  noimal  stress  resultants  are  given  by  the 


following  two  equations. 


cos  ijr 


I 

m=l 


As,  Nq  =  N.  +  (1+cos 
m  m  og  fg 


4)  \  A  S 

^  '  mm 


/ 

m=l 


(40) 


The  bending  moments  are  found  with  equations  (20)  using  the  following^ wal 1 
curvature  changes. 


4 


4 


V  = 


(1+cos  i) 
R 

o 


('-iD 


Again,  the  preceding  functions,  constants  Ajjj  and  A^,  and  computation  of  series 
coefficients  are  discussed  further  in  the  appendices. 

(b)  Particular  Solutions  -  Particular  solutions  applicable  to  the  nosecap 
depend  on  temperature  functions  (35)  and  (38)  that  vary  with  cos  \(i  in  the 
examples  of  Part  VIII.  These  solutions  are  written  below  in  closed  form. 
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A  cos 
o  ^  ’ 


=“A  'cos 
o 


♦ 


(42) 


They  are  equal  in  pairs,  and  A^  and  Ap  are  constants  that  are  evaluated  from 
temperature  data  in  Part  VI  and  in  the  appendices. 


4.  DISPLACEMENTS 

(a)  Wall  Extension  -  For  the  nosecap  of  Figure  2,  equation  (27)  reduces  to 


w  =  w^sin  i|(  -  (l+v)R^Q^/Eh  (43) 

where  is  defined  by  integral  (8)  and  computed  in  Reference  [27]. 

(b)  Radial  Displacement  -  The  radial  component  of  displacement  (u)  is 
obtained  from  equations  (28)  and  (43)  as 


u 


R 

o 

Eh 


(N  -  vN^  )  +  (1  +  v)N 


♦J 


w  cos 
0 


(44) 


where  the  coordinates  are  shown  in  Figure  2,  and  the  other  terms  are  given  in 
the  list  of  nomenclature.  The  foregoing  relations  for  u  and  lo  are  also 
developed  in  Appendix  G. 


VIII.  NUMERICAL  EXAMPLES 
1.  BICENTRIC-OGIVE  RADOME 

(a)  Thermal  Stresses  -  Based  on  stress  equations  presented  in  the  preceding 
text  and  the  appendices,  a  computer  program  was  set  up  for  stress  analysis  of 
the  bicentric-ogive  radome  described  in  Reference  [23]. 

The  computed  stresses  were  obtained  for  the  temperature  distribution  given 
by  equations (29)  and  (30)  using  the  outer-surface  straight-line  spanwise 
temperature  distribution  (T^)  shown  in  Figure  3  along  the  mainbody  as  written 
be  low . 


T  =  589  +  705  cos  i 
a  ^ 

Along  the  nosecap,  the  temperatures  were  represented  by 

T  =  50  +  1269  cos  A 
a  ^ 


(45) 


(46) 
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which  approximates  the  points  plotted  near  the  left  side  of  Figure  3  between 
z  =  28"  and  32". 

Computerized  thermal  stresses  are  listed  in  Tables  1  and  2  and  plotted 
on  Figure  6,  which  follow  immediately. 

(b)  Displacements  -  The  thermoelastic  displacements  were  computed  with 
jquations  (27)  and  (28)  for  the  mainbody  and  (43)  and  (44)  for  the  nosecap. 
The  numerical  results  are  tabulated  in  Table  3  and  illustrated  in  Figure  7. 

2.  RADOME  BASE  CONNECTIONS 


Four  kinds  of  attachment  of  a  radome  base  are  symbolically  illustrated 
in  Figure  8.  A  stiff  connection  is  indicated  by  8a  while  8b  is  flexible. 
The  attachments  represented  by  8c  and  8d  are  intermediate  between  the  stiff 
and  flexible  connections  diagrammed  in  8a  and  8b. 

3.  THICK-WALLED  HEMISPHERICAL  RADCME 


For  purposes  of  comparison,  a  thick-walled  (Ro/h  =  1.056)  hemispherical 
radome  was  analyzed  using  the  temperature  functions  given  by  equations  (29) 
to  (32),  inclusive.  Stress  values  from  the  computer  program  are  reported 
in  Tables  4  and  5  and  plotted  on  Figure  9.  Components  of  displacement  are 
listed  in  Table  6  and  shown  on  Figure  10, 

4.  THIN-WALLED  HEMISPHERICAL  RADOME 


In  the  third  example,  stresses  were  calculated  for  the  thin-walled 
(R  /h  =  26.5)  hemispherical  radome  sketched  on  Figure  11.  For  comparison 
wi?h  examples  one  and  two,  the  same  temperature  distributions  were  used. 
They  are  defined  by  equations  (29)  to  (32),  inclusive. 

The  thermal-stress  values  are  presented  in  Tables  7  and  8  and  shown 
on  Figure  11.  The  comi  ited  thermoelastic  displacements  are  listed  in  Table 
9  and  plotted  on  Figure  12. 


IX.  DISCUSSION 


1 .  NUMERICAL  EXAMPLES 

(a)  First  Example  -  The  thermal  stresses  and  displacements  that  are 
presented  in  Tables  1,  2,  and  3  for  the  bicentric-ogive  radome  are  illustrated 
in  Figures  6  and  7.  They  were  computed  with  v  =  0.244  and  E  =  16,400,000 
psi  as  mean  values  for  Pyroceram  9606  as  described  in  Appendix  A-4  and 
Reference  [l9].  Also,  the  straight-line  spanwise  temperature  distributions 
given  by  equations  (45)  and  (46)  were  employed  in  the  computations.  Thickness 
distributions  found  in  Reference  [22]  and  depicted  in  Figure  4  are  defined 
by  equations  (29)  and  (30). 
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Table  1.  Meridian  Stress,  0^,  psi  (Figure 


Coordinate  Ansle.  t.  Decrees 

1 

30 

64'’19'23"  68°3r56"  73°44'23" 

81°53'25" 

90 

+1.00  -33,580 

-32,910 

-23,302  -29, -425 

-28,615 

-25,121 

-29,448 

+0.75  -16,927 

•16,769 

-11,601  -18,304 

-18,124 

-16,123 

-19,995 

+  .50  -  6,470 

-  6,503 

-  2,944  -  9,885 

-  9,903 

-  8,695 

-11,398 

+  .25  +  3,699 

+  3,493 

+  5,621  -  1,552 

-  1,762 

-  1,338 

-  2,860 

0  +11,459 

+11,104 

+12,038  +  4,639 

+  4,245 

+  3,898 

+  3,571 

-  .25  +16,704 

+16,299 

+16,944  +  9,401 

+  8,924 

+  7,971 

+  9,009 

-  .50  +19,010 

+18,634 

+19,713  +12,075 

+  11,570 

+10,098 

+12,593 

-  .75  +18,457 

+18,239 

+20,997  +13,334 

+12,900 

+11,075 

+15,193 

+15.465 

+21.306  +13,676 

+13,380 

+11,305 

+17,157 

Table  2. 

Hoop  Stress,  Og,  psi  (Figure  6) 

_ Coordinate  Angle.  *.  Degrees _ 

c  0 

30 

64‘’19’23"  68°3r56"  73“44'23" 

8r53'25" 

90 

+1.00  -33,580 

-32,899 

-53,818  -30,254 

-28,901 

-25,445 

-31,846 

+0.75  -16,927 

-16,784 

-42,117  -19,022 

-18,494 

-16,397 

-23,715 

+  .50  -  6,470 

-  6,548 

-33,460  -10,905 

-10,356 

-  8,918 

-16,441 

+  .35  +  3,699 

+  3,411 

-24,895  -  2,042 

-  2,298 

-  1,509 

-  9,227 

0  +11,459 

+10,976 

-18,479  +  4,268 

+  3,627 

+  3,779 

+  4,122 

-  .25  +16,704 

+16,113 

-13,572  +  9,153 

+  8,225 

+  7,905 

11 

-  .50  +19,010 

+18,372 

-10,798  +11,953 

+10,790 

+10,086 

+  2,245 

-  9.210 _ +13.818 _ 

+12.441 

KSflREH 

+  4.149 

Table  3.  Thermoalastic  Displacements 

(Fiaure 

Z) _ 

(Deg-min-sec) 

cA 

(in)  (in  x  1(?)  (in  x  1C?  )  (In  x  ICp  ) 

(in  X  10?) 

0 

0 

28,175  23.6411  0 

0 

23.6411 

30 

0.132 

28,060  20.4585  -11.6288 

+0.1584 

23.5319 

64-19-23 

0.238 

28.025  9.8718  -20.9448 

-0.1785 

23.1540 

66-20-21 

1.188 

25.957  9.4295  -19.3686 

+0.8638 

21.5247 

68-31-56 

2.138 

23.671  8.6435  -17.8434 

1.5136 

19.7688 

70-57-49 

3.088 

21.096  7.8290  -16.2209 

2.1101 

17.8874 

73-44-23 

4.038 

18.110  6.8369  -14.4297 

2.5231 

15.7669 

75-55-48 

4.685 

15.724  6.3032  -13.0757 

2.9353 

14.2158 

78-31-17 

5.331 

12.871  5.5926  -11.5193 

3.1884 

12.4018 

81-53-25 

5.978 

9.124  4.7401  -  9.5873 

3.3403 

10.1601 

90 

6.625 

0  0  0 

0 

0 
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\ 

I 


Table  4. 

Meridian  Stress, 

,  psi  (Figure  9), 

1 

Coordinate  Anele.  ♦.  Deerees 

c 

0 

15 

30 

45 

60 

75 

90 

+1.00 

-33,878 

-34,194 

-35,6)5 

■  -40,348' 

-45,757- 

-  53;84y 

+0.75 

-18,333 

-18,775 

-20,166 

-22,687 

-26,600 

-32,216 

-  39,890 

+  .50 

-  8,954 

-  9,369 

-10,633 

-12,804 

-15,958 

-20,177 

-  25,591 

+  .25 

58 

-  406 

-  1,421 

-  3,017 

-  5,053 

-  7,349 

-  9,776 

0 

+  6,165 

+  5,936 

+  5,336 

+  4,619 

+  4,187 

+  4,527 

+  6,068 

-  .25 

+  9,542 

+  9,540 

+  9,692 

+10,452 

+12,524 

+16,753 

+  23,917 

-  .50 

+  9,454 

+  9,803 

+11,091 

+14,026 

+19,692 

+29,382 

+  44,302 

-  .75 

+  5,697 

+  6,591 

+  9,632 

+15,855 

+26,844 

+44,479 

+  70,516 

-1.00 

-  1,983 

-  257 

+  5,446 

+16,630 

+35,588 

+65,032 

+107,500 

Table  5, 

.  Hoop  Stress,  00,  psi  (Figure  9). 

X 

Coordinate  Angle.  4i.  Degrees 

c 

0 

15 

30 

45 

60 

75 

90 

+1.60 

-33,6'78 

-34,24'8" 

-35,324 

-TF,'93'7  '■ 

-38,732 

-39T()l!2  ■ 

'■■-F'/^TIUU 

+0.75 

-18,333 

-18,942 

-20,705 

-23,394 

-26,483 

-28,893 

-28,569 

+  .50 

-  8,954 

-  9,615 

-11,516 

-14,382 

-17,625 

-20,133 

-19,904 

+  .25 

58 

749 

-  2,723 

-  5,650 

-  8,882 

-11,277 

-10,915 

0 

+  6,165 

+  5,473 

+  3,516 

+  683 

-  2,315 

-  4,315 

-  3,513 

-  .25 

+  9,542 

+  8,924 

+  7,216 

+  4,860 

+  2,630 

+  1,679 

+  3,625 

-  .50 

+  9,454 

+  8,987 

+  7,755 

+  6,267 

+  5,354 

+  6,143 

+  9,982 

-  .75 

+  5,697 

+  5,503 

+  5,120 

+  5,132 

+  6,432 

+10,076 

+17,014 

-1.00 

-  1,983 

-  1,741 

770 

+  1,611 

+  6,369 

+14,447 

+26,012 

Table  6.  Thermoelastlc  Displacements  (Figure  10)- 


♦ 

(Deg .) 

u 

(In.  X  10*) 

w 

(in.  X  1(J*  ) 

“r 

(in.  X  lO*  ) 

(in.  X  id*  ) 

0 

4.798  . 

0 

0 

4.798 

15 

4.540 

-0.431 

0.759 

4.496 

30 

3.806 

-  .762 

1.243 

3.677 

45 

2.722 

-  .907 

1.283 

2.566 

60 

1.502 

-  .814 

0.894 

1.456 

75 

0.457 

-  .483 

0.317 

0.585 

90 

0 

0 

0 

0 

25 
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Table  7 

.  Meridian 

Stress,  I 

Psi  (Figure  11). 

1 

Coui'dinate  Angle. 

♦.  Decrees 

c 

0 

15 

. JQ. 

45 

60 

75 

+1.00 

-41,856 

-41,197 

-39,246 

-36,281 

-32,959 

-24,851 

-6,219 

+0.75 

-24,882 

-24,538 

-23,515 

-21,997 

-20,422 

-15,375 

-2,514 

+  .50 

-13,809 

-13,608 

-13,011 

-12,136 

-11,263 

-  8,069 

+  295 

+  .25 

-  2,869 

-  2,809 

-  2,630 

-  2,388 

-  2,200 

-  848 

+3,006 

0 

+  5,870 

+  5,793 

+  5,564 

+  5,188 

+  4,711 

+  4,235 

+2,569 

-  .25 

+12,426 

+12,251 

+11,732 

4-10,928 

■^10,035 

+  3,009 

+3, OS 9 

-  .50 

+16,489 

+16,239 

+15,494 

-^14,362 

+13,186 

+  9,75? 

+  733 

-  .75 

+18,388 

+18,103 

+17,249 

+15,978 

+14,765 

+10,207 

-2,657 

-1.00 

+18,886 

+18,591 

+17,705 

+16,416 

+15,326 

+  9,816 

-6,727 

Table 

8.  Hoop  S 

tress, 

psl  (Figure 

11). 

Coordinate  Angle.  ‘ 

(i.  Degrees 

c 

0 

15 

45 

60 

75 

90 

+1.00 

-41,856 

-41,197 

-39,270 

-36,308 

-32,246 

-25,066 

-28,906 

+0.75 

-24,882 

-24,537 

-23,535 

-22,035 

-19,830 

-15,077 

-22,150 

+  .50 

-13,809 

-13,607 

-13,025 

-12,186 

-10,792 

-  7,253 

-16,341 

+  .25 

-  2,869 

-  2,807 

-  2,640 

-  2,449 

-  1,850 

+  491 

-10,562 

0 

+  5,870 

+  5,795 

+  5,558 

+  5,115 

+  4,938 

+  6,102 

-  6,901 

-  .25 

+12,426 

+12,254 

+11,731 

+10,844 

+10,137 

+10,409 

-  4,245 

-  .50 

+16,489 

+16,242 

+15,498 

+14,266 

+13,163 

+12,694 

-  3,454 

-  .75 

+16,489 

+16,242 

+15,498 

+14,266 

+13,163 

+12,694 

-  3,659 

-1.00 

+18,886 

+18,596 

+17,719 

+16,296 

+15,049 

+13,845 

-  4,511 

Table  9. 

Thermoelastic  Dlspli 

icements  (Figure  12). 

u 

w 

Ur 

(in.  X  10" 

')  (In 

.  X  1(f) 

(In.  X  10* 

(in. 

X  1(f) 

0 

.  7.694 

0 

0 

7 

.694 

15 

+7.539 

-0. 

107 

+1.941 

7 

.285 

30 

+7.078 

-  . 

209 

+3.521 

6 

.140 

45 

+6.314 

-  . 

245 

■14.447 

4 

.482 

60 

+5.485 

-  . 

278 

+4.736 

2 

.767 

75 

+5.105 

-1. 

714 

+4.887 

1 

.487 

90 

-1.123 

0 

-1.123 

0 
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In  all  of  the  tables  and  figures,  negative  values  denote  compressive 
stresses  and  positive  values  are  tensile  stresses.  Component  u  (Figure  8) 
of  displacement  is  positive  radially  outward  along  the  R-axis,  and  w  is 
positive  with  increasing  iji.  For  the  rectangular  components  of  displace¬ 
ment,  w^  is  positive  in  the  positive  direction  of  the  Z-axis,  and  u^  is 
positive  in  the  direction  of  increasing  r  which  rotates  about  the  Z-axis 
as  located  by  the  coordinate  angle  8. 

As  shown  on  Figures  6  and  7  the  bicentric-ogive  radome  is  fixed  at 
the  base.  Various  kinds  of  radome-missile  attachments  with  regard  to 
joint  stiffness  or  flexibility  are  illustrated  in  symbolic  form  on  Figure 

8. 


From  the  numerical  results,  it  is  seen  that  the  stresses  generally 
decrease  from  the  tip  to  the  base  of  the  radome.  The  maximum  tensile 
hoop  stress  of  19,010  psi  at  the  tip  of  the  nosecap  is  less  than  the 
allowable  material  strength  of  35,000  psi.  In  the  mainbody,  the  meridian 
stress  along  the  inner  surface  of  the  wall  decreased  from  +21,310  psi  at 
64“19'23"  to  +13,680  psi  at  68°31'56".  This  finding  is  consistent  with 
the  18,000  psi  calculated  from  measured  strains  between  these  two  interior 
points  as  referred  to  in  Part  VI-2b  on  strength  data. 

(b)  Second  Example  -  For  the  thick-walled  hemispherical  radome  in 
this  example,  the  base  attachment  is  shown  by  Detail  "a"  of  Figure  8, 
and  the  dimensions  are  given  in  Figures  9  and  10. 

From  the  numerical  results  in  Tables  4  and  5,  it  is  seen  that  the 
stresses  increase  from  the  tip  to  the  base  of  the  radome  while  the  converse 
occurred  in  the  previous  example.  These  stresses  are  excessive  and  would 
crack  a  real  radome.  Distributions  of  these  thermal  stresses  are  plotted 
in  Figure  9  for  comparison  only.  And  the  displaced  position  of  the  radome 
wall  caused  by  aerodynamic  heating  is  sketched  in  Figure  10. 

(c)  Third  Example  -  The  base  connection  is  represented  by  Detail  "b" 
of  Figure  8  for  the  thin-walled  hemispherical  radome  in  this  example,  and 
the  dimensions  are  shown  in  Figures  11  and  12.  Also,  the  radome  wall  has 
the  same  thickness  and  its  central-surface  radius  is  the  same  as  the  base 
radius  of  the  bicentric-ogive  radome  in  the  first  example. 

From  the  tabulated  numerical  results  (Tables  7  and  8  and  Figure  11), 
one  observes  that  the  maximum  stresses  occur  at  the  tip.  There  is  18,890 
psi  tension  along  the  inner  surface  and  41,860  psi  compression  along  the 
outer  surface  of  the  wall  at  the  tip  of  the  radome.  Meridian  tensile 
stress  maxima  decrease  from  tip  to  base  along  the  span  and  fall  inside 
the  wall  beyond  45“,  while  the  compressive  maxima  remain  along  the  wall's 
outer  surface  and  decrease  along  the  span  to  6220  psi  at  the  base. 

Again,  as  shown  here  in  Figure  12,  the  wall  expands  outward  over  most 
of  the  length,  but  is  induced  by  the  spanwise  heating  to  move  inward  at 
the  base.  With  the  stiff  base  attachments  employed  in  examples  one 
(Figure  7)  and  two  (Figure  10),  such  inward  movement  is  suppressed. 
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2 .  SUMMARY 

(a)  General  Remarks  -  With  fixed  ends,  examples  one  and  two,  displace" 
ment  and  wall-slope  change  are  prevented  at  the  base  of  a  radome ;  i.e.,  the 
boundary  conditions  imposed  by  the  radome-missile  joint  are:  Vjj  =  uj,  =  wjj  = 

0,  which  represent  conditions  that  are  more  severe  than  occur  in  actual 
structures.  A  missile  body  expands  and  contracts  with  changing  temperatures, 
so  ub  would  not  vanish . 

The  bicentric-ogive  radome  wall  material  was  found  to  sustain  the  thermal 
stresses  by  both  test  and  theory.  Of  the  two  hemispherical  radomes  that 
were  analyzed  for  purposes  of  comparison,  the  thick-walled  radome  is  over¬ 
stressed,  but  stresses  in  the  thin-walled  hemisphere  lie  within  safe  limits. 

(b)  Recommendations  -  It  appears  that  the  methods  of  analysis  of 
thermal  stresses  and  displacements  in  ogival  radomes  under  aerodynamic 
heating  that  are  presented  herein  are  satisfactory  for  numerical  evaluations. 
Further  study  of  the  problem  is  needed,  however,  to  provide  answers  to  un¬ 
resolved  questions. 

It  is  therefore  recommended  that  the  investigations  be  continued  to 
obtain  information  on  nose  temperatures,  stress  maxima,  and  the  effect  of 
different  kinds  of  radome-missile  connections. 


X.  CONCLUSIONS 


Based  on  the  reported  investigations,  it  is  concluded  that: 

(1)  The  methods  of  analysis  are  satisfactory  for  the  ca’  jlation  of  thermal 
stress  and  displacement  in  ogival  radomes. 

(2)  Additional  information  is  required  on  temperature  distributions. 

(3)  Further  studies  are  necessary  on  nosecap  behavior  and  thermal  stress 
maxima . 

(4)  Further  investigation  of  radome-missile  joint  effects  is  desirable. 


APPENDIX A:  THEORY 


1.  EQUILIBRIUM  REQUIREMENTS 

(a)  Stress  Equations  of  Equilibrium  -  The  volume  element  in  Figure  1 
is  oriented  with  respect  to  meridian  lines  (located  from  the  center  of 
curvature  by  radius  R)  and  the  angle  ijr,  which  measures  the  inclination  of 
R.  The  third  coordinate  is  9,  which  is  the  rotational  angle  about  the 
geometric  Z-axis  for  the  ogival  body  of  revolution.  They  are  often  referred 
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to  as  the  natural  ogive  coordinates.  Orthogonal  coordinate  systems  (R,  9,  \|i) 
and  (r,  6,  z)  are  shown  in  Figure  la. 

For  axially  symmetric  bodies  and  heat  distributions,  the  partial 
differential  equations  of  equilibrium  for  a  volume  element  follow  immediately. 


2i(rRa  )  a(rT  ,) 


dy 


+  =  ra^  +  ROgsin 


a(rR^T  )  a(rRCT,) 


ay 


Jl 


a^ii 


\1;  . 
■*—  =  R  CpCos 


(47) 


Equations  (47)  are  derived  from  E  F  =  0  and  E  F^  =  0, 
geometrical  and  thermal  symmetry,  tnere  are  no  changes 


and  T  .  and  t  -  are  zero  as  well  as  E  F 


yG 


\)i9 


e 


0. 


respectively.  Owing  to 
in  the  0  direction 


ILL.  Stress  and  Moment  Resultants  -  When  equations  (47)  are  integrated 
over  the  wall  thickness,  the  equations  of  equilibrium  are  obtained  in  terms 
of  the  stress  and  moment  resultants  of  integrals  (6),  (8),  and  (10). 


dt|i 


(r  N,  +  R  NflSin  )|f)-(rRp  +rRp 
'  -  CO  '  a  a  av  "  " 


C  'ij 


s  s  sy 


(48) 


d(r  N  ) 

— =  (R  N-cos  \|i-rQ,)-(rRp,+rRp,) 
d\ji  c  9  ^  c  \|(  a  a  ai|(  s  s  s\|i 


(49) 


d(r  M  ) 

—  ^  =  R  (M-cos  \(r  +  rQ)  +  c(rRp  -  rRp  )  (50) 

d\l(  c  9  ^  c^^  s  s  sijt  a  a  a\|i 

In  the  foregoing  relations,  p^y  and  p  ,  are  the  radial  and  tangential  components 
of  pressure  at  the  outer  surface  of  the  radome;  and  Pgy  and  Pg^,  at  the  inner 
surface.  These  pressures  are  equated  to  zero  in  a  thermal -stress  problem. 

For  the  investigation  of  temperature  stresses,  therefore,  the  resultant 
stress  and  moment  conditions  that  ensure  equilibrium  are  expressed  as  follows. 

d(r  Q. /R  sin  iji)  d(r  M  ) 

N^=Q^coti,  Ng  =  - =  Mgcos  ^  +  r^Q^  (51) 

These  expressions  hold  for  both  thick  and  thin-walled  ogive  shells  from  the 
inner  apex  to  attachment  points  at  the  missile  body,  but  not  within  the 
attachment  region  owing  to  restraints  set  up  by  the  connections.  In  any  case. 
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however,  the  body-joint  problem  requires  further  analysis  and  can  often  be 
treated  by  reintroducing  the  surface  pressures  as  derived  in  equations  (48) 
to  (50),  inclusive. 

2,  STRAIN  EXPRESSIONS 

(a)  Thermal  Strain  -  For  Pyroceram  9606,  the  thermal-strain  function 
(Cj.)  was  developed  from  test  data  in  Reference  [l9].  From  these  results 
together  with  temperature  distributions  reported  in  Reference  [22]  for 
wind  tunnel  testing  of  a  Von  Karman  radome,  the  function  is  defined  as 
follows . 


=  (-241.2  +  3.445T)  x  lO' 


70°  ^  T  <  400°F 


=  (-119.1  +  2.958T)  X  10-®,  200°  s  T  ^  600°F 


=  (+190,3  +  2.221T)  X  lO'® ,  300°  ^  T  ^  1400°F 


Temperature  (T)  is  given  by  equations  (29)  to  (31),  inclusive. 

Inspection  of  the  distribution  of  temperature  defined  throughout  the  radome 
by  these  equations  and  Figure  3  revealed  that  equation  (52)  satisfactorily 
describes  the  temperature  in  the  inside  half  of  the  wall  (-1  ^  x  ^  0)  along 
the  entire  span;  equation  (53)  defines  T  over  the  full  span  in  the  quarter 
thickness  above  the  wall  centerline  (0  s  x  s  +0.5);  and  (54)  serves  the  same 
purpose  in  the  outer  quarter  of  the  wall  thickness  (+0.5  ^  x  +1).  This 
arrangement  expedites  the  calculation  of  temperature  functions,  integrals  (7) 
and  (11),  and  the  stress  equation  (2). 

The  strains  computed  with  equation  (52)  agree  with  the  measured  data. 
Reference  [l9],  within  -1  to  +3  percent.  Those  calculated  from  equation  (53) 
correspond  to  the  measurements  within  five  percent;  and  the  ones  found  with 
(54)  agree  within  seven  percent  in  the  low  part  of  the  temperature  range,  which 
is  seldom  used,  and  deviations  in  the  upper  half  of  the  range  lie  between 
approximately  -2  and  +4  percent. 


(b)  Actual  Strain  -  Components  of  strain  for  elastic  behavior  were 
thoroughly  discussed  by  Love  [l]  and  other  authors  [2],  [3],  [4],  and  [5], 
We  have 


'e  '  <'ec  ■ 

1  dw  1  / 

e  =;r~  n  +  —  ,  =  —  u  sin  il’  +  w  cos  ili 

i|[c  R^  dill  J  ec  r^  (  ^ 
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V 


du 


d^  ’ 


V  cos  ill 

^0  “  r 

c 


(57) 


where  equations  (55)  define  the  normal-strain  components  at  any  point  in  the 
wall  of  the  radome ,  equations  (56)  give  these  strains  at  the  central  surface, 
and  (57)  shows  the  change  (V)  of  the  slope  of  the  wall  together  with  curvature 
changes.  General  solutions  for  the  latter  functions  appear  as  equations  (22); 
and  particular  solutions,  as  equations  (25)  and  (26). 

(c)  Elastic  Strain  -  As  the  result  of  thermal  strain  due  to  temperature 
change  plus  elastic  strain  produced  by  constraints  against  thermal  expansion 
or  contraction,  we  have 


e 


'1> 


't’ 


(58) 


where  e^,  eg  are  the  components  of  elastic  strain  in  the  wall  of  the  heated 
radome.  And  the  other  terms  are  given  by  equations  (55)  to  (57),  inclusive. 

3.  MECHANICAL  PROPERTIES 

(a)  Young's  Modulus  -  Dependence  on  temperature  of  the  mechanical 
properties  of  Pyroceram  9606  was  examined  in  Reference  [l9].  In  the  pre¬ 
sent  numerical  examples,  an  average  value  of  E  was  employed  and  is  written 
below. 


E  =  E  =  16,400,000  psi;  70°  S  T  ^  1400°F 
av  ’  ’ 


(59) 


Values  of  E  over  the  above  temperature  range,  which  includes  all  temperatures 
in  the  examples,  agree  within  four  percent  of  the  quoted  average. 

(b)  Poisson's  Ratio  -  As  pointed  out  in  Reference  [l9],  Poisson's  ratio 
(v)  is  nearly  constant  for  the  radome-wall  material.  For  the  specified 
temperature  range,  its  mean  value  is 


V  =  V  =  0,244,  70°  ^  T  S  1400°F  (60) 

av  ’ 

and  the  largest  differences  between  tnis  mean  value  and  the  test  data  amount 
to  ±2.5  percent. 

4.  STRESS  STRAIN  RELATIONS 

r 

(a)  Hooke's  Law  -  The  previously  discussed  wall  stresses  and  strains 
are  related  by  means  of  Hooke's  law;  e.g., 


i5 
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Ee  =0  -  va-,  Eej,  =  a»  -  va,  (6; 

\Jf  t|/  6  0  6  iji 

from  which  th“  stresses  are  customarily  expressed  in  the  following  manner. 


“e  ■  + 'e>-  e'  -  e/(i  -  v=)  (6: 

As  presented  in  the  next  two  articles,  the  preceding  relationships  let  us 
express  the  central-surface  strains  and  curvature  changes  in  terms  of  the 
stress  and  moment  resultants. 

(b)  Normal-Stress  Resultants  -  When  E  and  \)  are  represented  by  their 
constant  mean  values,  equations  (59)  and  (60)  respectively,  the  integrals 
(5),  (6),  and  (7)  together  with  equations  (55)  to  (58)  and  (62)  lead  to 
the  following  expressions  for  and  Ng . 


Second  and  higher  order  effects  are  neglected  in  this  calculation  of  and 
Ng  in  terms  of  the  central-surface  strain  components. 

(c)  Bending  Moments  -  Equations  (9),  (10),  and  (11)  are  used  with  (55) 
to  (61),  inclusive,  to  obtain 


Ml-  -D(x  +  vXq),  1^=-D(vx  +Xe) 


wherein  the  second-order  effects  are  again  neglected,  D  is  the  wall  flexural 
stiffness  factor  stated  by  equation  (21),  and  the  curvature-change  components 
are  given  by  equations  (22)  and  (57). 


APPENDIX  B;  TEMPERATURE  FUNCTIONS 


1.  NORMAL-STRESS  TEMPERATURE  FUNCTICTIS 

Functions  that  depend  on  the  temperature  distribution  and  have  the 
dimensions  of  a  ncrmal-stress  resultant  are  defined  by  equations  (7)  and 
can  be  expressed  in  the  forms  written  below. 


Nte  =  %+  KiCOS 


(65) 
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When  equations  (2),  (29),  (30),  (31),  (52),  (53),  (54),  (59),  and  (60)  are 
substituted  into  (7)  and  integrated,  the  are  found  to  have  the  following 
values . 


Ko=  4.648Eh/(l-v)  x  lO^  , 

iq=  6.792Eh/(l-v)  X  l(f 

(66) 

Kg-  4.393Eh/(l-v)  x  lO'  , 

1^=  3.420Eh/(l-v)  X  lO’^ 

(67) 

The  radii  rj,  and  R^.  that  occur  in  equation  (65)  are  illustrated  in  Figure  1 
and  defined  by  equations  shown  on  the  figure.  They  are  evaluated  in  Reference 
[25]. 

2.  BENDING -MOMENT  TEMPERATURE  FUNCTIONS 


Functions  that  depend  on  the  temperature,  distribution  and  have  the 
dimensions  of  a  bending -moment  resultant  are  defined  by  equations  (11)  and 
can  be  expressed  in  the  forms  written  below. 


\e  ^  '^0+  +  (Jg+  J3C0S  (68) 

When  equations  (2),  (29),  (30),  (31),  (52),  (53),  (54),  (5C , ,  and  (60)  are 
substituted  into  (11)  and  integrated,  the  are  found  to  have  the  following 
values  , 


Jo=  1.269Eh^/(l-v)  X  10*, 

Ji=  0.989Eh^/(l-v)  X  10* 

(69) 

Jj,=  1.704Eh^ /(1-v)  X  10% 

J^=  2.195Eh^/U-v)  X  iO"^ 

(70) 

One  can  see  that  the  have  the  dimensions  of  pounds  per  inch  (ppi)  ; 
J.,  of  inch-pounds  per  inch  (ippi).  Again  the  above  are  evaluated  in 

[is]. 


and  the 
Reference 


APPENDIX  C :  DIFFERENTIAL  EQUATIONS  FOR  SHEAR  AND  SLOPE  FUNCTIONS 


1.  AUXILIARY  SHEAR  FUNCTION 


To  find  solutions  for  the  total  shear  (Q^),  equation  (8),  the  auxiliary 
shear  function  (q)  is  introduced  such  that 


j7 
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R  q  cos  ijt 

q  =  r  Q,/R  sin  \|r,  N  =  - 

^  c  i|i  c  ^  ’  \|i 

and,  by  solving  equations  (63),  the  central-surface  normal  strains  are  obtained 
in  terms  of  q. 

From  equations  (56)  and  (57),  V  is  found  in  terms  of  the  central-surface 
normal  strains. 


N  --a 
diji 


(71) 


CSC  \li  *^^^c^9c^ 

Rc  d^ 


e  cot  ill 

\j|C  ^ 


(7?.) 


And  with  equations  (5),  (63),  and  (71)  put  into  (72),  the  following  ordinary 
differential  equation  is  obtained 


(L^+  v)q  =  EhV  +  (N  .  -  vN  .)cot  ^  ^ 


'till  “te^ 


R  dr^ 


vN  ) 

tilt 


(73) 


The  operator  is  defined  as  follows: 

r  CSC  ij)  2 


L^q  = 


R  df 


+  cot  l|l 


d\|(  '1' 


(74) 


and,  in  general,  the  solutions  of  equation  (73)  are  obtained  as  infinite 
series  for  q  and  V,  which  are  exemplified  by  equations  (12)  and  obtained 
in  Reference  [26]. 


2.  WALL- SLOPE  CHANGE 

Relationships  similar  to  the  foregoing  are  derived  for  V  by  substitution 
of  equations  (11),  (22),  and  (23)  into  (51).  The  result  is 


(L^-  v)(EhV)  =  (B  +  v^) 


M  -cot  ill  ,  d(r  M  ) 

to  ^  CSC  i!;  c  t'lt 

R  "  "  R3  diji 

c  c  ^ 


(75) 


where  is  given  by  equation  (74),  q  by  (71),  and  fhe  temperature  moments 
by  (11)  and  (68),  The  wall-bending  parameter  (B)  Is  calculated  with 
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which  is  developed  from  the  expression  EhR^/D,  where  D  is  given  by  equation 
(21).  The  general  solutions  of  equation  (75)  are  given  in  Reference  [26]. 

3.  SPECIFIC  RESULTS 

When  specific  temperature  functions  defined  by  equations  (65)  and  (68) 
are  put  into  equations  (73)  and  (75),  the  results  are  as  follows. 


(L^+  v)(q  +  1^)  =  EhV  +  (l-v)iqr^/R^  +  KgR^/r^tan  i)t  (77) 


(L^-  v)(EhV)  =  (B*+  V^)(Jir^/R^  +  *4 /^c  " 

Particular  solutions  of  (77)  and  (78)  are  necessitated  by  the  presence  of 
constants  K^,  and  general  solutions  when  these  constants  are  absent. 
These  solutions  are  developed  in  Reference  [24]  . 


where  series  f^  (m  =  1,  2,  3,  4)  are  given  by  equations  (12)  and  (16). 

The  series  coefficients  (3,^^^  computed  from  recurrence  formulas  as 

explained  in  Reference  [26],  and  constants  B  and  B'  are  related  as  follows. 

m  m 
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Bi=  (vB^+  3^Bj,)/Eh,  Bg'=  (vBjg-  3"Bi)/Eh  (81) 

And  similar  relations  hold  for  m  =  3  and  4. 

(b)  Particular  Solutions  -  With  reference  to  equations  (77)  and  (78), 
only  the  terms  involving  r^,  as  the  only  variable  can  be  solved  in  closed  form. 
These  forms  are  asssociated  with  Kj  and  .  We  have 


(L=+  v)qp  =  EhVp  +  (l-v)Kir^/R^  (82) 

0^-  v)(EhV  )  =  (8*+  v^)(J^r  /R^  -  q  )  ^  (83) 

r  L  L.  L' 


and  the  particular  solutions  are  written  below. 


=  B-r  /R  , 
6  c  C 


(84) 


Jj  (l-v2)(Ji+  R^K^) 

^  “  F  ■  a+F)r 


>  I  _ 


(l-v)(v®+  8*)(Ji+  R  K^) 


(l+p4 )R^Eh 


(85) 


These  solutions,  equations  (84),  are  included  in  equations  (23)  and  (25). 
For  other  terms  (the  ones  that  involve  ,  and  Jg  ) ,  however.  Infinite 

expansions  are  required  to  satisfy  the  differential  equations.  In  the 
current  examples,  approximation  functions  are  used  that  meet  the  require¬ 
ments  with  errors  of  one-half  percent  or  less.  The  approximate  solutions 
are  discussed  in  Reference  [24]. 


APPENDiy  jundary  conditions 


U _ GENERAL  REQUIREMENTS 


In  general,  one  considers  the  requirements  at  the  juncture  of  two 
arbitrary  segments  of  the  mainbody,  such  as  MBl  and  MB2  with  the  junction 
at  \|i-  =  73°44'23"  in  the  first  numerical  example.  Eight  functions  are 
involved:  displacements  u,  w  and  slope  change  V;  stress  resultants  Q^j,, 
N^,  N^ ;  and  moments  M^,  Mq.  At  any  juncture  (ijf^)  where  there  is  wall 
material  and  geometric  continuity,  all  eight  functional  requirements  can 
be  met  by  the  following  five  equations. 
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w(f^+)  =  w(il/^-l). 

v(i^+)  =  v(^.-) 

(86) 

q('l/^+)  =  q('lt^-i)» 

Nq()^^+)  =  Nq(i^^-) 

(87) 

(88) 

From  these,  one  constant  in  w(\lf£+)  is  determined  directly  and  four  relations 
for  constants  of  integration  are  obtained.  The  equations  and  evaluation  of 
the  constants  are  discussed  in  Reference  [26]. 

2.  RADOME  BASE 


At  the  base  of  a  radome,  three  equations  are  found  according  to  the  type 
of  attachment  as  illustrated  in  Figure  8.  In  all  cases,  axial  relative  dis¬ 
placement  between  missile  and  radome  vanishes,  or 

w(90^)  =  Wj^  =  0  (89) 

and  the  constant  w^  in  Wj^  is  computed  directly.  Then  for  a  fixed  base,  we 
have 


u(90°)  =  Uj^  =  0,  V(90“)  =  =  0  (90) 

and  for  a  free  base,  the  requirements  are  stated  below. 


q(90°)  =  =  0, 


(91) 


With  a  hinged  base,  one  uses 


u 


b 


=  0, 


(92) 


and  with  a  bearing-supported  base 


q 


b 


=  0, 


(93) 


where  the  relationship  between  q  and  Q 


is  given  in  equation  (71). 
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3.  NOSECAP-MA INBODY  JUNCTION 


At  =  64 “19 '23",  which  locates  the  interface  between  the  nosecap  and 
MBl  (Figure  2)  in  the  bicentric-ogive  radome  of  the  first  example,  the 
requirements  for  strain  continuity  lead  to  the  following  expressions. 


(94) 

Xq(^-i+)  =  Xq(’\-), 

»dX^('^i+)  =  RgX^('!-i-) 

(95) 

('I'l  -)  + 

eQc^'^i-)  +  ^Xe(^i-)  =  0 

(96) 

Equations  (94)  and  (96)  produce  satisfaction  of 


N3(^,+)  =  N2(^,-)  (97) 

which  also  meet  similar  requirements  imposed  on  and  u,  since  the  constant 
Wg  in  the  nosecap  formula  for  w  is  computed  with  the  next  equation. 


w('lti+)  =  w(\;^  -)  (98) 

The  first  of  equations  (95)  satisfies  the  requirement  on  V,  equation  (86), 
and  the  second  one  fulfills  a  requirement  that  is  similar  to  equation  (88). 


APPENDIX  F;  THERMOELASTIC  DISPLACEMENTS 


1.  WALL  EXTENSION 


From  equations  (56)  for  the  mainbody,  we  have 
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))rc 


r  €  A/  N.  -  vNg  r  (vN.  -  H, ) 

c  9c  ^  sin  d(w  esc  \lt)  _  ^  c  ^ 

R  sin  ill  R  d\li  Eh  R  Eh  sin  iji 

c  ^ 


d\|t 


R  Eh  sin 
c  ^ 


Eh 


r  (vN  , 
c  t\li 


R  Eh  sin  ^ 


^  ,  (r  +  vR  sin  A  )  , 

+  fR  sin  J,  +  vr  )  ^  ^ - £ - £ - —  ^ 

(R^s  ij(  r^Eh  sin  di)? 


(99) 


which  is  integrated  to  secure  equation  (27).  The  integration  constant  in 
equation  (27)  is  w^,  and  functions  w^  and  are  shown  below. 


w  =  W  sin  ill 

q 


(r^  +  vR^sin  i|i)q 
Eh  sin  \(f 


(100) 


w 


t 


dill 

sin  ill 


(101) 


W 


q(2r^ 


+  R^sin  i|i)cos  I  dill 
r^  siit^iji 


(102) 


When  the  normal- Jtress  temperature  functions  in  equation  (101)  are  given  by 
equations  (65),  w^  becomes 


'^t  "  "^h^  cos  ill  +  (Kgcos  ill  +  ) 

R  sin  ill 


+ 


Eh  sin 


^  1'^ 


{\+  Ajsec  ^  )  +1^/^ 


Aj  =  -Ln  (esc  ill  -  cot  ill), 


1  +  cos  ( ilH-i|i  ) 

- - —  ,  /^=  'tn(R^/r^csc  1^) 

c  c  ; 


(103) 

(104) 


where  the  \  are  logarithmic  expressions  as  indicated  above.  The  trigonometric 
sine  of  the  central-surface  apex  is  computed  with 


sin  ill  =  H/R  ,  H  =  R  -  r, 
^c  c  c  be 


(105) 
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and  sec  \(i^  can  then  be  readily  obtained.  Values  of  the  Kj  are  calculated 
from  equations  (66)  and  (67).  Functions  W,  integral  (102),  can  be  computed 
with  infinite  series,  but  it  is  usually  more  convenient  to  find  W  by  finite 
summation. 

2.  RADIAL  DISPLACEMENT 


To  find  u  for  the  mainbody,  wc  again  use  the  second  of  equations  (56) 
together  with  equation  (27)  for  w.  When  is  obtained  from  equation  (63),  the 
following  result  is  found. 


u  =  (Ng  -  vNj^  )r^/Eh  sin  ijr  -  w  cot  \[t  (106) 


Equation  (106)  can  be  used  throughout  the  ogival  mainbody  since  \(t  >  0. 


APPENDIX  G:  NOSECAP  EXPRESSIONS 


1 .  TEMPERATURE  FUNCTIONS 

(a)  Normal-Stress  Temperature  Functions  -  When  one  uses  outer-surface 
temperature  distributions  that  vary  with  cos  \(i  in  the  manner  of  equations  (45) 
and  (46)  for  the  nosecap,  the  normal-stress-resultant  temperature  function, 
equation  (35),  can  be  expressed  as  follows. 


"t  *  "tt  ' 

For  the  nosecap  on  the  bicentric-ogive  radome ,  surface  temperatures  were 
given  by  equation  (46),  and 


I^=  -0.0612Eh/(l-v)  X  Id*,  13.873Eh/(l-v)  x  lO*  (1 

which  were  employed  in  the  first  example. 

For  the  hemispherical  radomes,  the  surface  temperatures  were  given  by 
equation  (45),  and 


Kp'=  5.842Eh/(l-v)  X  lO*  ,  1^=  7.772Eh/(l-v)  x  l(f  (109) 

which  were  used  in  the  second  and  third  examples  . 
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(b)  Bending-Moment  Temperature  Functions  -  Since  temperature  functions 
for  the  bending -moment  resultant  have  the  same  general  form  as  equation  (107), 
one  can  use 


“e  ■  '  ”te  ■  •'d  +  •’i' 

and,  for  the  nosecap  with  surface  temperatures  from  equation  (46),  the  con¬ 
stants  are  written  below. 


Jq'  =  0.5195Eh^/(l-v)  X  lO*  ,  =  2 .830Eh® /(1-v)  x  l(f  (111) 


For  the  hemispherical  radomes,  the  surface  temperatures  were  given  by 
equation  (45),  and 


Jg'  =  1.732Eh^/(l-v)  X  10*,  =  1.585Eh^/(l-v)  x  10* 


(112) 


which  were  employed  in  the  second  and  third  numerical  examples. 

2 .  DIFFERENTIAL  EQUATIONS  FOR  SHEAR  AND  SLOPE  FUNCTIONS 

(a)  Total  Shear  -  For  hemispherical  segments,  such  as  the  nosecap,  equation 
(73)  reduces  to 


(l2+  v)Q^  =  EhV  -  (1-v)  ^  (113) 

and  equation  (74)  becomes 

'''  d^  "  ^ 

which  lead  to  general  solutions  of  the  kind  indicated  by  equations  (39)  and 
(40)  that  are  more  fully  described  in  Reference  [27]. 

(b)  Wall-Slope  Change  -  For  the  nosecap,  equation  (75)  reduces  to 

r  ,  dM  ~ 

tt"-  v)(Ehv)  -  -<f^  +  ^)=)  +:^-^  (115) 

L  -J 

where  L"  is  given  by  equation  (114),  by  (8)  and  the  temperature  moment  by 
(38)  and  (110).  The  general  solutions  are  discussed  in  Reference  [27]. 
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(c)  Specific  Results  -  When  the  specific  temperature  functions  defined 
by  equations  (107)  and  (110)  are  put  into  equations  (113)  and  (115),  the 
results  are  as  follows. 


(L^+  v)  =  EhV  +  (l-v)K^'  sin  i)t  (116) 

(L^-  v)(EhV)  =  (^+  v^) 

Particular  solutions  of  these  equations  are  necessitated  by  the  presence  of 
constants  K^'  and  and  general  solutions  when  these  constants  are  absent. 

3.  SOLUTION  OF  EQUATIONS 

(a)  General  Solutions  -  To  find  general  solutions  of  equations  (116)  and 
(117),  one  observes  that 


Jj'sin  \|( 


-  Q,i 


(117) 


a"+  Bo)  Qg  =  0, 


(L*+  Vg  =  0 


(118) 


where  subscript  "g"  denotes  that  the  solutions  are  general.  They  are  found  to 
be 

4  4 


Q  =  sin  il( 


A  a  , 
m  m 


V  =  sin  Jt 
g  ^ 


A'a 
m  m 


m=l 


m=l 


(119) 


where  the  infinite  series  (s  )  are  indicated  by  the  first  of  equations  (39), 
and  they  are  developed  in  Reference  [27]. 

(b)  Particular  Solutions  -  The  terms  associated  with  sin  in  equations 
(116)  and  (117)  possess  closed-form  solutions.  They  are 


Qp  =  Apsin  ii, 


=  -A^sin  I 


(120) 


where  constants  ^  ,  Ag  depend  on  K^'  ,  as  shewn  below. 

J^'  (l-v=)(j;+  R^K^')  (l-v)(v^+  0‘^)(J^'+  R^K^) 

(1  +  B„^)R  Eh 


(1  +  B  )R 


A'  = 
o 


(121) 


Numerical  values  of  the  foregoing  constants  were  obtained  with  the  aid  of  equations 
(108),  (109),  (111),  and  (112)  and  used  in  the  three  example  problems. 
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4.  BOUNDARY  CONDITIONS 

(a)  NosecaP  Requirements  -  The  nosecap  boundary  conditions  at  its 
juncture  with  the  mainbody  are  discussed  in  Appendix  E-3  in  conjunction  with 
equations  (94)  to  (98),  inclusive.  However,  owing  to  errors  in  the  original 
computer  program  not  all  of  these  conditions  were  met  in  the  first  numerical 
example.  The  numerical  results  reported  herein  satisfy  requiirri'-nts  that 
meridian  fiber  stresses  at  the  inner  and  outer  surfaces  of  the  wall  are 
equal  between  mainbody  and  nosecap,  but  there  are  discontinuities  in  some 

of  the  terms  that  appear  in  equations  (94)  to  (97),  inclusive.  Equation  (98) 
was  exactly  fulfilled  however. 

In  consequence,  the  nosecap  stresses  must  be  regarded  as  approximate 
rather  than  precise  results.  Furthermore,  they  are  believed  to  be  reasonably 
good  approximations. 

(b)  Base  Connection  Requirements  -  The  base-connection  boundary  conditions 
given  by  equations  (89)  to  (93)  were  used  in  comparison  studies  with  hemi¬ 
spherical  radomes.  They  were  completely  satisfied  in  these  studies.  In 

the  second  example  in  this  report,  conditions  (90)  were  imposed;  and  in 
the  third  examnle,  conditions  (91). 

5.  DISPL/iC 

(a)  as  ion  -  Extension  (w)  of  the  nosecap  wall,  which  is 

Illustrated  i  g^re  8,  is  computed  from  equation  (56)  with  r^  =  sin  ; 
i  .c . , 


yc 


-  e 


Gc 


sin 

R 


CSC  ,^)^  (1,±_V2  /M  .  M  ^ 

d^  Eh  ^  1  2^ 


.  (L+Jil  sin  4  diSLc^c^jlL) 

Eh  ^  dill 


(122) 


.rron  which  the  following  expression  is  obtained. 

w  =  w^  sin  ^  -  (l+v)R^Q/Eh  (123) 

And  w  is  an  arbitrary  constant, 
o  •' 

(b)  Radial  Displacement  -  Again  using  the  central-surface  hoop  strain  and 
equations  (56)  and  (63),  one  has 


^0^  =  (u  +  w  cot  ’1i)/Rq  =  (Ng-  vN^  )/Eh 


'  1 


(124) 
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which  is  solved  to  obtain  the  following  expression  for  u. 


u  =  (Nj-  vNj)R^/Eh  -  V  cot  ^  (125) 

Since  cot  \ji  becomes  infinite  at  =  0,  it  is  desirable  to  substitute  equation 
(123)  into  (125)  for  numerical  work.  In  this  way,  we  arrive  at  the  next 
relationship  for  the  radial  displacement  (u). 


u 


R 

o 

Eh 


(Ng-  vNi  )  +  (1  +  v) 


w  cos  \li 
o  ^ 


(126) 


The  preceding  equation  for  u  was  reported  earlier,  equation  (44),  as  being  the 
most  suitable  one  for  the  computer  program. 
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The  solution  that  is  developed  herein  can  be  used  to  compute  thermal 
stresses  and  displacements  in  the  moderately  thick  walls  of  appreciably  curved 
ogival  shells  whose  material  properties  change  with  heat  intensity  along  the  span 
and  through  the  wall.  No  restriction  is  placed  on  temperature  dependent  vari¬ 
ations  of  the  thermal  stress  expression. 

Three  examples  of  the  computerized  results  are  presented.  The  first 
example  deals  with  a  28.  3"  x  6.  75"  bicentric-ogive  radome;  the  second,  with  a 
thick-walled  hemispheric  radome;  and  the  third,  with  a  thin-walled  hemisphere. 

In  these  calculations,  linear  spanwise  temperature  distributions  are  used  that 
were  obtained  from  experimental  and  heat-transfer  analyses  of  a  radome  tested 
in  a  wind  tunnel.  Through  the  waD,  over  specified  layers,  the  thermal- strain 
function  is  approximated  as  being  linear  with  temperature.  However,  tiie  tern 
perature  varies  as  a  third- order  nonlinear  function  throu^  the  thickness  of  the 
wall  as  found  for  the  test  radome.  It  appears  that  errors  due  to  the  approximation 
amount  to  less  than  four  percent  in  general. 
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